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Abstract 

We establish a connection between the general equations of nonlin- 
ear elastodynamics and the nonlinear ordinary differential equation of 
Pinney [Proc. Amer. Math. Soc. 1 (1950) 681]. As a starting point, 
we use the exact travelling wave solutions of nonlinear elasticity dis- 
covered by Carroll [Acta Mechanica 3 (1967) 167]. The connection 
provides a method for finding new exact and approximate dynamic 
solutions for neo-Hookean and Mooney-Rivlin solids, and for the gen- 
eral third- and fourth-order elasticity models of incompressible solids. 
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1 Introduction 



The equations of elastodynamics belong, under the usual constitutive as- 
sumptions, to the class of hyperbolic systems. It is known that if a hyperbolic 
system is nonlinear, then smooth solutions of initial- value or initial boundary- 
value problems do not usually exist globally in time, and that singularities 
will develop, typically after a finite time, even when the initial or bound- 
ary data are smooth pQ. Existence of smooth solutions to the initial- value 
problem of nonlinear elastodynamics is possible only in special situations and 
only for body waves. 

To the best of our knowledge, the only explicit examples of smooth so- 
lutions to the equations of nonlinear elastodynamics are the Carroll waves 
[21 El SI IS]- These solutions were first investigated in [2] in the form of finite- 
amplitude, circularly-polarized, shear plane waves traveling in any equi- 
biaxially deformed solid, as 

x = p,sX + A cos(kZ — ojt), y = figY + A sin(kZ — ut), z = \$Z, (1.1) 

where A, k, u, us and Xs are arbitrary constants (in incompressible materials, 



Recently these solutions have been generalized by Destrade and Sacco- 
mandi [HI Ej to dissipative and dispersive materials and in the case of a time 
dependent homogeneous underlying motion by Rajagopal [8] and Pucci and 
Saccomandi [9]. 

The aim of this note is to show that there exists a nice connection between 
some of the remarkable solutions of Carroll and the Milne-Pinney equation 
ordinary differential equation (also known as the Ermakov equation [TU]). 
This is the nonlinear ordinary differential equation 



for c constant and p(x) given. In such a way it is possible to uncover some sim- 
ple new exact solutions to the equations of nonlinear elastodynamics. These 
new solutions show once again the relevance of the Carroll waves solutions 
to our understanding of elastodynamics. 

2 Basic Equations 

We call x (X, t) the current position of a particle which was located at X in 
the reference configuration. We follow Pucci and Saccomandi [9] and consider 
the motion 



x = ii{t)X + f(Z,t), y = n{t)Y + g(Z,t), z = X(t)Z, (2.1) 



/4A S = i). 



y" + p{x)y + cy 3 = o, 
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which consists of two shearing motions / and g in the Z-direction, combined 
with a time- dependent homogeneous biaxial stretch A(t) (here, p = A -1 / 2 ). 
This class of motions has been considered by Rajagopal [8J; it is a general- 
ization of the Carroll waves (1.1). 

Two kinematic quantities associated with this motion are the deformation 
gradient and the left Cauchy-Green strain tensor, 

F = dx/dX, B = FF T , (2.2) 

respectively. We focus on incompressible solids, which can undergo only 
isochoric motions, so that det F = 1 at all times (this is indeed the case for 
this motion). 

We also require that the material be hyperelastic and isotropic, and so 
we introduce the strain energy density W = W(Ii,I 2 ), where 1\ and l 2 are 
the first and second principal invariants of B, respectively. For isochoric 
motions, they are given by 

h = tr B, I 2 = tr(B^). (2.3) 

The general representation formula for the Cauchy stress tensor T is 

T = -pi + 2W 1 B - 2W 2 B-\ (2.4) 

where p is the yet indeterminate Lagrange multiplier introduced by the con- 
straint of incompressibility, and W\ = dW/dli, W 2 = dW/dI 2 . 

Now, the balance equation of linear momentum, in the absence of body 
forces, is 

div T = p d 2 x/dt 2 , (2.5) 

where p is the (constant) mass density. 

Standard computations [9J lead to the following specialization of the bal- 
ance equations to the motion ( 2.1[ ), 

p (f tt + p tt X) = -p x + [2{XW 1 + W 2 )f z \ z , 

P {gu + IhtY) = -p y + \ [2(AWi + W 2 )g z } z , 

p\ tt Z = -p z + [2A 2 ^! - 2A- 1 (/i 2 /l + g 2 z ) W 2 ] z , (2.6) 

where letter subscripts denote partial differentiation. Then we write these 
equations in the reference configuration, and we introduce the shear strains 
F = fz and G = gz, and the generalized shear modulus Q = 2(W\ + A -1 !^)- 
We point out that here 

h = 2p 2 + A 2 + F 2 + G 2 , I 2 = 2\ 2 + \-\p 2 + F 2 + G 2 ). (2.7) 
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Finally we eliminate the Lagrange multiplier p, and we obtain the deter- 
mining equations for the strains F and G, 



P\F U 



Ptt 



{QF) 



zz ' 



P 



G t 



Ptt 



G 



(QG) 



zz ' 



(2.1 



where Q = Q(F 2 + G 2 ,t). We assume that the Baker- Ericksen inequalities 
hold, so that Q > at all times 

It is convenient to recast (2.8) as a single differential equation by intro- 
ducing the complex quantity A = F + iG, to obtain 



P 



fi 



[QM 



ZZ ' 



where now Q = Q(A 2 ,t). 



We seek solutions of (2.9) in the form 

A = fi(t) exp[i(kZ + <f>(t))] 



(2.9) 



(2.10) 



where the wave-number k is a real constant, and Q (the amplitude) and </> 
(the phase) are unknown real functions of time. Notice that this type of 
solutions gives 



f(z,t) 



sin (kZ + <f)(t)), g(Z,t) 



~k~ 



cos (kZ + 4>(t)) , 

(2.11) 

for the shear motions, making the connection with the Carroll waves (1.1): 
here we have the same spatial variations, and potentially richer temporal 
variations. 

Now the full set of partial equations (2.9) reduces to the following system 
of ordinary differential equations, 



y + £ 



- k 



2 Q 



n. 



b"Q + 20! c 



0. 



(JL p 

in the unknowns Q(t) and (f)(t) (now Q 
system in the unstretched case, when \i 



(2.12) 



Q(n 2 ,t)). Carroll [3] derived this 
A = 1. 

= h/n 2 , 



Direct integration of the second equation in (2.12) gives 0' 
where k± is an arbitrary constant, which we take to be positive without 
loss of generality. The amplitude Q(t) is now the solution of the following 
nonlinear and non-autonomous ordinary differential equation, 



k 



(2.13) 



ji p 

The case k± = (giving 0(t) = constant) has been discussed in detail by 
Pucci and Saccomandi [9], who establish a clear relationship with Melde's 
problem of parametric resonance. 
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3 The Milne-Pinney Equation 



We first specialize the equations to the Mooney-Rivlin solid, with strain- 
energy density 

W(h, I 2 ) = |// (| + (3) (h - 3) + lfi {\ - 0) (J 2 - 3), (3.1) 

where /io > is the infinitesimal shear modulus and /3 is a constant such 
that —1/2 < (3 < 1/2 (when /3 = 1/2 we recover the special case of the 
neo-Hookean material). In that case, the generalized shear modulus is 

Q(t) = ^ [I + / 3+(I- / 3)A- 1 (t)], (3.2) 



which is independent of Q, and the governing equation (2.13) becomes 



n „_j^_^ [|+/3 + (i _ /3)A - 1] j fi _| = 



(3.3) 



This is indeed the Milne-Pinney or Ermakov equation. A detailed commen- 
tary about this equation is given in a recent paper by Leach and Andriopolous 

HDJ. 

Note that this is not the first time that this differential equation arises 
in the framework of nonlinear elastodynamics. In fact, it has arisen pre- 
viously in the study of large amplitude oscillations of thin-walled tubes of 
Mooney-Rivlin materials. This problem was first treated by Knowles [12] for 
a general hyperelastic incompressible material, and the connection between 
this problem and the Milne-Pinney equation was made by Shahinpoor and 
Nowinski [13], see also the book by Rogers and Ames [13] . Notice however 
that there is a noteworthy difference between the results in [13] and the 
results here: the oscillations of a Mooney-Rivlin tube are governed by the 
Milne-Pinney equation in the approximation of a thin-walled tube, whereas 
here no approximation at all was necessary to derive the equation. This fact 



was alluded to by Carroll in [3] who noticed the reduction (2.11), but did 
not use it to deduce some exact solution. 

To fix ideas, take the following initial conditions 

ft(0) = Vl > 0, ft'(0) = 0. (3.4) 



Then the general solution of (3.3) can be obtained by using two linearly 



independent solutions u and v, say, of the linear equation 



n" - - ^ [| + p + (| - 0) a- 1 ] } n = o, (3.5) 
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such that 

w(0) = tt , u'(0) = 0, 



and 



v (o) = o, v '(o) = i/n - 



(3.6) 



As shown by Pinney |15j . the general solution of (3.3) is then 



Q= ^u 2 + k\v 2 . (3.7) 
An example can be worked out when \(t) is a constant, static equi-biaxial 



pre-stretch, X = \ s , say. Then (3.3) reduces to 

2 _ f^ok 2 ri 



a" + ujIvl + 



0, where w n 



[i + /3+(i-/3)A s 1 ] >0. (3.8) 



tt 3 p 
In this case, u = cos(ojot), v = [sm(coot)} /(^o^o), an d the Pinney solution 



(3.7) is 



k 2 

fi(t) = \ I 'fig cos 2 (u Q t) + sm 2 (oj t). 
From the equation <fr' = ki/Q 2 we deduce the phase to be 



(f)(t) = tan 



-i 



k\ , \ 
— -2 tan(w t) 



^0- 



(3.9) 



(3.10) 



where 0o is the initial phase. 

More complex solutions may be obtained by applying the approach of 
Shahinpoor and Nowinski [13J to equation ( 3.3[ ). For simplicity of exposi- 
tion, we consider the neo-Hookean case, /3 = 1/2. Then, the Milne-Pinney 



equation (3.3) reduces to 

Q" - 



p" pok 2 
p p 



k 2 



(3.11) 



Now consider the case where the neo-Hookean solid is stretched linearly with 
time, \(t) = 1 + at, say, where a > is a constant. Then p" / p = 3a 2 /[4(l + 
at) 2 }. The following changes of function and of variable, 



co 



and 



C = a(l + at), where a = 




give a non-dimensional version of (3.11) 



3 \ 

Utf + 1 - ~ ~ I CO 



4C 2 



where 



c = 



co c 



Po k 

p a ' 



p k 2 tt^ 



(3.12) 



(3.13) 
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with initial conditions changed from (3.4) to 

cu(a) = 1, 



<J {a) = 0. 



(3.14) 



The linear Pinney equation corresponding to (3.13) is a Bessel equation, and 



we express the set of fundamental solutions in terms of Bessel functions of 
the first and second kind as 



u, v 



(vr/4) [AjJt (C) + Bfr (C)] , j = 1, 2, 



(3.15) 



where 



A 1 = -Y 1 (a) + 2aY \ 
B\ = Ji(a) — 2«J (q 



a) 



A 2 = -2aYi(a), 
B 2 = 2aJ x (a). 



(3.16) 



These constants have been computed in order to accommodate the initial 
conditions for u and v, which are now: 



u(a) 



u'(at) = 0, 



and 



v(ce) = 0, v'(a) = 1. 



(3.17) 



The resulting non-dimensional solution is u 



\u 2 + cv 2 ] 1 ^ 2 . Its behavior 



depends on two constants: c > defined in (3.13), which is arbitrary because 



k\ is arbitrary; and a > defined in (3.12), which is equal to the ratio of the 



speed of infinitesimal shear waves to the speed of stretching. This connection 
leads us to define a 'supersonic' range, < a < 1 and a 'subsonic' range, 
a > 1. 

At the origin ( = a, we have (3.14) and also u>^(a) = c — 1 + 3/(4a 2 ). 
Hence, when a is in the supersonic range < a < ^/3/A ~ 0.866, the 
amplitude u always grows at first because ou^(a) = 0, uj^(a) > 0. Otherwise, 
when a is in the supersonic range a/3/4 < a < 1 or in the subsonic range 
a > 1, the amplitude can either grow or decay, depending on the value of 
c with respect to 1 — 3/(4a 2 ). In Figure [l] we show the early variations of 
uj when c = 1.0, a = 0.8 (thick plot), c = 1.0, a = 2.0 (medium thickness 
curve), and c = 0.5, a = 2.0 (thin curve). In all cases, the amplitude 
oscillates and remains bounded over time. 



4 The Full Non-Linear Equation 



We established that the Milne-Pinney equation emerges from the specializa- 



tion of the full non-linear equation (2.13) to the Mooney-Rivlin strain energy 



density (3.1). In fact, the exact solutions obtained subsequently by the Pin- 



ney method are valid for the whole class of general third-order non-linear 
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Figure V. Variations with time of the amplitude of a generalized Carroll wave 
in a neo-Hookean solid under a linearly growing stretch (non-dimensional 
units). 

incompressible solids. These are described by the following expansion of the 
strain energy W, 

W = no tr (E 2 ) + ^tr (E 3 ) , (4.1) 

where E = (F T F — I)/2 is the Green strain tensor and A is a third-order 
elasticity constant in the notation of Hamilton et al. [16J. Indeed, it is 
known since Rivlin and Saunders [17J that when the Mooney-Rivlin strain 
energy (3.1) is expressed in terms of tr (-E 2 ) and tr(i£ 3 ), and neglecting 



higher powers, it embraces (4.1), once the following identification is made: 
A = 2/io(2/3 — 3). It follows that a solution valid for Mooney-Rivlin solids is 
also valid for all third-order incompressible solids (but the reverse is not true 
in general). 

Further, we may use the exact Pinney solutions for Mooney-Rivlin solids 
as the starting point of a perturbation scheme in the fourth-order elasticity 
theory of incompressible solids. As shown by Ogden |l8j . only three elastic 
constants are necessary to describe these. In the notation of Hamilton et al. 
6], the expansion of the strain energy density reads 

W = Hotr (E 2 ) + 4 tr (E 3 ) + D (tr(E 2 )) 2 , (4.2) 
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where A and D are third-, and fourth-order elasticity constants, respectively. 
From Destrade et al. [21] we deduce that the following strain energy 



W = ii Q (§ + /?) (l! - 3) + (|-/?)(/ 2 



3) + ^(/i-3) 2 



(4.3) 



covers (|4.2|), once terms of fifth-order in E and higher are neglected, and the 

A/2 + D 



following identifications apply. 





3 A 
2 + W 



7 = 1 + 



/'o 



(4.4) 



Using (4.3), the corresponding generalized shear modulus Q is easily com- 
puted as 



Q = Po [\ + P + [\ - 0) A" 1 + 7(2/i 2 + A 2 + n 2 - 3)] . 



(4.5) 



As an example, consider a small-amplitude vibration of period u, super- 
imposed upon a large static stretch X s , so that A(t) = As[l + ecos(u;£)], where 
|e| << 1. Then we expand (2.13) up to term of order e 1 . We find 

Hok 2 



--r(2Xs 1 + a| - 3 - n 2 )n 



p 

, Pok 2 
H e cos 



( U t)[(±-(3)\ s 1 + 2 1 (\; s 1 -\ 2 s )]n, (4.6) 



where c^o is defined in (3.8). 

It is possible to consider several perturbation schemes that exploit the 
solution (3.9) as the zeroth-order approximation of equation (4.6). The sim- 
plest of such perturbation scheme is obtained considering that the last term 
in (4.3) is a perturbation of the Mooney-Rivlin strain energy, i.e. that 7 = 76, 
where 7 = 0(1). With that assumption, we find that with 

k\ /x fc 2 



n 3 



p 



n" + u 2 Q- 7 ± = [(\-0) Xs 1 cos M) - li 2 ^ 1 + A| - 3 - Q 2 )} Q. 

(4.7) 

We may then perform a standard multiple scale expansion [19] starting from 
the exact solution (3.9) of the Milne-Pinney equation at e = 0. In this way 
it is possible to investigate parametric resonance phenomena such as those 
revealed by the Melde string experiment [20]. This has been done in the 
linear case by Lord Rayleigh [22] , and more recently in a complete nonlinear 
context by Pucci and Saccomandi [9] (the main difference is that here we 
have traveling waves instead of standing waves.) 
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5 Concluding Remarks 



The discovery of Carroll waves has had a long-lasting influence in the field of 
wave propagation in nonlinear elasticity. These solutions are elegant exam- 
ples of exact reductions, and are also explicit examples of exact, closed-form 
solutions of an highly non-linear theory of solid mechanics. Moreover, the 
peculiarities of these solutions proved fundamental to a better understand- 
ing of the physical and mathematical properties of the classical theory of 
elastodynamics. Fritz John |23j was the first to recognize this fact. 

To give an idea of the recent research originating from the papers by 
Carroll [21 S IH [5] , we mention: the note by Rubin and Rosenau [21] in the 
framework of the theory of a nonlinear string; the papers by Destrade and 
Saccomandi [5J 125] on rotating media, viscoelastic materials of differen- 
tial type, and dispersive non-linear elastic materials, respectively; the paper 
by Rajagopal [B]; and the paper by Pucci and Saccomandi [9] on Melde's 
phenomena. 

The aim of this note was to point out a link between the Carroll waves 
and the Milne-Pinney equation. It is clear that the solutions presented are 
somewhat puzzling from a mathematical point of view. Indeed, we know 
that global existence of body waves is guaranteed in nonlinear elasticity if 
the null condition is satisfied [26], but this result holds only for small initial 
data. The Carroll wave (1.1), and its various generalizations, is an explicit 
example of global existence of solutions for any value of the amplitude of 
initial data. It is still necessary the investigate the mathematical structure 
lying hidden behind such simple solutions. 



Acknowledgments 

The work of MD is supported by a Marie Curie Fellowship for Career Devel- 
opment awarded by the Seventh Framework Programme of European Com- 
mission. The work of GS is partially supported by Italian GNFM. 



References 

[1] Dafermos, C, Hyperbolic Conservation Laws in Continuum Physics, 
Second Edition, Springer, Heidelberg, 2005. 

[2] Carroll, M.M., Some results on finite amplitude elastic waves, Acta Me- 
chanica 3 (1967) 167-181. 



10 



[3] Carroll, M.M., Oscillatory shearing of nonlinearly elastic solids, ZAMP 
25, (1974) 83-88. 

[4] Carroll, M.M., Plane circular shearing of incompressible fluids and 
solids, Q. J. Mech. Appl. Math. 30 (1976) 223-234. 

[5] Carroll, M.M., Finite amplitude standing waves in compressible elastic 
solids, J. Elasticity 8 (1978) 323-328. 

[6] Destrade, M. and Saccomandi, C, On finite amplitude elastic waves 
propagating in compressible solids, Phys. Rev. E, 72 (2005) #0016620. 

[7] Destrade, M. and Saccomandi, C, Nonlinear transverse waves in de- 
formed dispersive solids, Wave Motion 45 (2008) 325-336. 

[8] Rajagopal, K.R., On a class of elastodynamic motions in a neo-Hookean 
elastic solid, Int. J. Non-Linear Mech. 33 (1998) 397-405. 

[9] Pucci, E. and Saccomandi, C, Parametric resonance in non-linear elas- 
todynamic, Int. J. Non-Linear Mech. 44 (2009), 560-569. 

[10] Leach, P. G. L. and Andriapoulos, K., The Ermakov equation a commen- 
tary, Applicable Analysis and Discrete Mathematics 2 (2008) 146-157. 

[11] Beatty, M.F. Topics in finite elasticity: Hyperelasticity of rubber, elas- 
tomers, and biological tissues-with examples, Appl. Mech. Rev. 40 
(1987) 1699-1734. 

[12] Knowles, J. K., Large amplitude oscillations of a tube of incompressible 
elastic material, Quart. Appl. Math. 18 (1960) 71-77. 

[13] Shahinpoor, M. and Nowinski J.L, Exact solution to the problem of 
forced large amplitude radial oscillations of a thin hyperelastic tube, 
Int. J. Non-linear Mech. 6 (1971) 193-208. 

[14] Rogers, C and Ames, W.F, Nonlinear boundary values in science and 
engineering, Academic Press, 1989. 

[15] Pinney, E., The nonlinear differential equation y"(x) +p(x)y + cy~ 3 = 0, 
Proc. Amer. Math. Soc, 1950, 681. 

[16] Hamilton, M.F., Ilinskii, Y.A., Zabolotskaya, E.A., Separation of com- 
pressibility and shear deformation in the elastic energy density. J. 
Acoust. Soc. Am. 116 (2004) 41-44. 



11 



[17] Rivlin, R.S., Saunders, D.W., Large elastic deformations of isotropic 
materials. VII. Experiments on the deformation of rubber. Phil. Trans. 
Roy. Soc. A 243 (1951) 251-288. 

[18] Ogden, R.W., On isotropic tensors and elastic moduli. Proc. Cambr. 
Phil. Soc, 75 (1974) 427-436. 

[19] Stoker, J.J., Nonlinear vibrations in mechanical and electrical systems, 
Interscience, 1950. 

[20] Melde, F., Uber die Erregung stehender Wellen eines fadenformigen 
Korpers, Ann. Physik Chemie. 109 (1860) 193-215. 

[21] Destrade, M., Gilchrist, M.D., and Murphy, J.G., Onset of non-linearity 
in the elastic bending of beams, ASME J. Appl. Mech. (to appear 2010). 

[22] Strutt, J.W. (Lord Rayleigh), The Theory of Sound, 2nd. ed., Dover, 
1945. 

[23] John, F., Wave propagation in Trends in Applications of Pure Mathe- 
matics to Mechanics: v. 1 ([Monographs and studies in mathematics) 
edited by Gaetano Fichera, Pitman London (1976) 175-185. 

[24] Rosenau, P., Rubin, M. B., Some nonlinear three-dimensional motions 
of an elastic string. Phys. D 19 (1986), 433-439. 

[25] Destrade, M. and Saccomandi, G., Some results on finite amplitude 
elastic waves propagating in rotating media, Acta Mechanica, 173 (2004) 
19-31. 

[26] Domanski, W. and Ogden, R.W., On the null condition for nonlinearly 
elastic solids, Arch. Mech. 58 (2006) 339-361. 



12 



